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ABSTRACT 

Spitzer observations of extended dust in two optically normal elliptical galaxies provide a new confirmation of 
buoyant feedback outflow in the hot gas atmospheres around these galaxies. AGN feedback energy is required to 
prevent wholesale cooling and star formation in these group-centered galaxies. In NGC 5044 we observe interstellar 
(presumably PAH) emission at 8 fj,m out to about 5 kpc. Both NGC 5044 and 4636 have extended 70 /im emission 
from cold dust exceeding that expected from stellar mass loss. The sputtering lifetime of this extended dust in the 
~ IkeV interstellar gas, ~ lO*" yrs, establishes the time when the dust first entered the hot gas. Evidently the extended 
dust originated in dusty disks or clouds, commonly observed in elliptical galaxy cores, that were disrupted, heated 
and buoyantly transported outward. The surviving central dust in NGC 5044 and 4636 has been disrupted into many 
small filaments. It is remarkable that the asymmetrically extended 8 /im emission in NGC 5044 is spatially coincident 
with Halpha+[NII] emission from warm gas. A calculation shows that dust-assisted cooling in buoyant hot gas moving 
out from the galactic core can cool within a few kpc in about ^ 10^ yrs, explaining the optical line emission observed. 
The X-ray images of both galaxies are disturbed. All timescales for transient activity - restoration of equilibrium and 
buoyant transport in the hot gas, dynamics of surviving dust fragments, and dust sputtering - are consistent with a 
central release of feedback energy in both galaxies about lO'' yrs ago. 

Subject headings: galaxies: elliptical and lenticular; galaxies: ISM; infrared: galaxies; infrared: ISM 



1. INTRODUCTION 

In a recent survey of elliptical galaxies observed with 
the Spitzer telescope we found spatially extended cold in- 
terstellar dust emitting at 70 fim in a significant fraction 
of our sample (Temi, Brighenti & Mathews 2007). A few 
unusual galaxies in our sample, having atypically large 
and massive dust lanes that often dominate the optical 
image, probably acquired their dust and cold gas from a 
merger. However, some optically normal elliptical galax- 
ies also contain extended cold dust, such as the two galax- 
ies we discuss here, NGC 5044 and 4636. Evidently, the 
source of this dust is internal, not from mergers. Since 
this dust is thought to be in direct contact with the hot, 
virialized interstellar gas {kT ~ 1 keV), it has a short 
Hfetime ('^ 10'' yrs) to sputtering destruction by ther- 
mal ions. Consequently, this dust is a spatial tracer of 
extremely transient events that recently occurred on kpc 
scales. Both galaxies are extended at 70 /im, but NGC 
5044 is also extended at 8 /im in a manner similar to 
the highly asymmetric optical line emission from warm 
gas in this galaxy. The temperature of the warm gas is 
maintained by photoionization by UV radiation from hot 
post-asymptotic giant branch (AGB) stars and radiation 
of optical line emission; the dust is heated by starlight 
and hot thermal electrons and cooled by infrared emis- 
sion. 

As we proposed in Temi, Brighenti & Mathews (2007), 
current evidence is consistent with an internal origin for 
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this dust, and we hypothesize that dust has been buoy- 
antly transported from the galactic cores out to several 
kpc following a feedback heating event. In both galax- 
ies small, disorganized fragments of optically absorbing 
dusty gas are visible within the central '^100 parsecs, fur- 
ther evidence of a recent central energy release. In addi- 
tion, Chandra X-ray images reveal that the extended hot 
interstellar gas in both galaxies is undergoing short-term 
gasdynamical activity. All this information converges to 
confirm the transient state of affairs in these galaxies. 
In particular, the association of interstellar PAH emis- 
sion (from polycyclic aromatic hydrocarbon molecules) 
and warm gas (T - lO" K) in NGC 5044 indicates that 
we may be viewing this galaxy at a rare moment im- 
mediately following a release of energy near the central 
black hole. This type of central feedback energy release, 
thought to be caused by the accretion of a small mass of 
gas into the central black hole, has been proposed as a 
likely solution of the so-called cooling flow problem. In 
the following we adopt distances of 15 Mpc (1' = 4.36 
kpc) for NGC 4636 and 33 Mpc (1' = 9.60 kpc) for NGC 
5044 (Tonry et al. 2001). 

2. IMAGES OF NGC 5044 AND NGC 4636 

Figures 1 and 2 show Spitzer images of NGC 5044 
and NGC 4636 at 160, 70, 24, and 8 ^m. The 8-4.5 
fim images are differences of Spitzer images at 8.0 and 
4.5 fim. Details of the acquisition and reduction of far- 
infrared data at 24, 70 and 160 /im are described in Temi, 
Brighenti & Mathews (2007). These data were obtained 
with the Multiband Imager Photometer (MIPS) (Rieke 
et al. 2004) on the Spitzer Space Telescope. Both NGC 
5044 and 4636 are detected at 160 /im, but neither image 
can be distinguished from a point source. However, as 
we illustrate in Figure 10 of Temi, Brighenti & Mathews 
(2007), at 70 /um both NGC 5044 and 4636 are clearly 
extended beyond the MIPS point response function. 
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Mid-infrared data at 3.6, 4.5, 5.8 and 8.0 for NGC 
4636 were obtained using the Infrared Array Camera 
(IRAC) onboard Spitzer under the guaranteed time pro- 
gram (PID 69) led by PI G. Fazio (Fazio et al. 2004). 
NGC 4636 was observed for 60 seconds on each of the 
four IRAC wavelength channels. NGC 5044 was ob- 
served during an IRAC calibration program (PID 1160, 
PI: W. T. Reach). Images of NGC 5044 at each of the 
four IRAC bandpasses were recorded at 5 different po- 
sitions slightly offset from the galaxy center. Observa- 
tions at each position were acquired using a five posi- 
tion gaussian dither parameter with an integration time 
of 12 seconds per frame. Thus, the total on target in- 
tegration time at IRAC wavelengths was 300 seconds. 
For all wavelengths and with both MIPS and IRAC, full 
coverage imaging was obtained for all observations with 
additional sky coverage to properly evaluate the back- 
ground emission. We used the Basic Calibrated Data 
(BCD) products from the Spitzer Science pipeline (ver- 
sion 13.2) to construct mosaic images when necessary. 
Pipeline reduction and post-BCD processing using the 
MOPEX software package provide all necessary steps to 
process individual frames: dark subtraction, flat-fielding, 
mux-bleed correction, flux calibration, correction of focal 
plane geometrical distortion, and cosmic ray rejection. 

Also shown in Figures 1 and 2 are images of the 
IIa+[NII] emission from the warm gas in NGC 5044 
(isophotes) and 4636 (gray scale). The panels at the 
lower right in both figures illustrate isophotes of the X- 
ray images of the same central regions of both galaxies, 
produced from data in the Chandra archives, superim- 
posed on optical images from the Digital Sky Survey. 

The Ha+[NII] image of NGC 5044 is quite unusual, 
with extensions toward the East, North and (especially) 
to the South. These warm gas features are obviously 
unrelated to the underlying stellar distribution. The 
Ha+[NII] image in NGC 4636 is more compact, but nev- 
ertheless exhibits irregularities in the outer regions. The 
radial velocities of the warm gas in both galaxies de- 
termined from slit spectra show highly irregular activ- 
ity of ± ~ 150 km on scales of ~ 1 kpc that are 
unrelated to the bulk motion of tlw^ stars (Caon et al. 
2000). These quasi-random velocities are consistent with 
the disordered spatial appearance of the warm gas. 

3. EXCESS DUST IN NGC 5044 AND 4636 AND ITS ORIGIN 

Recent studies of mid-infrared Spitzer IRS spectra of 
elliptical galaxies indicate that most of the emission 
shortward of 8 /xm is photospheric, while emission in the 
8-20 fim region is from circumstellar dust around mass- 
losing red giant stars (Bregman, Temi & Bregman 2006; 
Bressan et al. 2006). In addition, Temi, Brighenti & 
Mathews (2007) find a tight correlation between fluxes at 
24 /zm from Spitzer MIPS data and optical B-band fluxes, 
suggesting that circumstellar dust also dominates at this 
wavelength. But MIPS fluxes from elliptical galaxies at 
70 and 160 /im are uncorrelated with optical fluxes from 
starlight, so we regard this colder dust as truly interstel- 
lar. It is remarkable that the fluxes at 70 and 160 /J,m 
vary by factors of 30-100 among galaxies with similar op- 
tical fluxes. Some elliptical galaxies with large 70 and 160 
/im fluxes contain very extended, multi-kpc disks of cold, 
dusty gas consistent with the strength and spatial extent 
of the 70 yum image. Evidently, these galaxies have expe- 



rienced an unusual major merger with a gas-rich spiral. 
However, a few optically normal elliptical galaxies, in- 
cluding NGC 5044 and 4636, also show unexpected large 
70 and 160 /um fluxes with 70 /um emission extended well 
beyond the MIPS point response function. 

Because the 8-24 ^m spectrum of elliptical galaxies 
is dominated by circumstellar dust in outflowing stellar 
winds, some internally produced dust is expected when 
this dust moves into the interstellar environment, where 
it is heated by diffuse starlight to much lower tempera- 
tures and radiates in the far-infrared. In Temi, Brighenti 
& Mathews (2007) we describe a simple steady state 
model in which interstellar dust is continuously produced 
by normal mass loss from a single population of evolv- 
ing stars and continuously destroyed by ion sputtering. 
A far-infrared spectral energy distribution (SED) can be 
calculated assuming a power-law initial grain size dis- 
tribution (X a~^'^ in < a < a^nax- This grain size 
distribution is normalized so that S = (l/150)z(r) is the 
ratio of dust to gas mass in the circumstellar outflows 
where z{r) « (r/i?e)~°'^'"' is a typical radial variation of 
the metal abundance (in solar units) in elliptical galax- 
ies (Arimoto et al. 1997). We assume amorphous silicate 
grains and compute the temperature of each dust grain. 
The grain production rate is proportional to the local 
stellar density and the specific stellar mass loss rate for 
an old stellar population a* = A'/*/A'/* = 4.7 x 10^^" 
s"-'^ (Mathews 1989). The sputtering rate of the grains is 
determined by the local density and temperature in the 
hot gas as determined from X-ray observations. The to- 
tal galactic SED is found by integrating the dust emission 
over the steady state dust size distribution and galactic 
radius. We show in Figure 11 of Temi, Brighenti & Math- 
ews (2007) that the far-infrared SED at 70 and 160 /im 
observed in at least a few galaxies - such as NGC 1399, 
1404 and 4472 - can be fully explained by this type of 
normal stellar mass loss. Elliptical galaxies that agree 
with our steady state dust model are among those with 
the lowest 70 and 160 /um luminosities. 

However, as illustrated in Figure 3, the 70 and 160 /xm 
fluxes of both NGC 5044 and 4636 are far in excess of 
our predicted steady state SED. (The 24 /im flux shown 
in Figure 3 is circumstellar and is unrelated to our SED 
model.) The mass and spatial extent of the "excess" in- 
terstellar dust in these galaxies can be estimated using 
our steady state model but artificially increasing the stel- 
lar mass loss rate within radius r^x- Increasing a* by 
some factor increases the mass of interstellar dust and 
rex can be adjusted until the mean grain temperature ap- 
proximately fits the 70 /Ltm/160 /im flux ratio observed 
in each galaxy. Interstellar dust is heated by starlight 
and by inelastic collisions with thermal electrons. Since 
the mean intensity of starlight and the electron density 
both decrease with galactic radius, the dust tempera- 
ture T4 also decreases, so rgx can be adjusted to change 
mean dust grain temperature (T^) and therefore the 70 
/im/160 /im flux ratio. These values of /„, rex and the 
excess dust mass estimated from them are only approx- 
imations, since it is unlikely that the radial distribution 
of the excess non-stellar dust exactly follows the stellar 
density profile, but with a sharp truncation at rex- 

Each panel of Figure 3 shows three combinations of 
fa and rex, holding Umax = 0.3/im constant. The best- 
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fitting SED for the extra dust in NGC 5044 and 4636, 
shown by dotted lines in both panels, corresponds respec- 
tively to a mass of excess dust of 1.5 x 10^ and 1.1 x 10^ 
Mq. The physical size of this dust, Tex ~ 4 — 5 kpc, 
is comparable to or exceeds the 70 ^m point response 
function for both galaxies. We propose below that r^x 
is limited by the sputtering time in buoyant dusty gas 
flowing from galaxy cores. 

We argue in Temi, Brighenti & Mathews (2007) that 
the incidence of (otherwise normal) elliptical galaxies 
with excess dust emission at 70/xm exceeds the fraction 
of ellipticals that are likely to have experienced mergers 
with dusty gas-rich (dwarf) galaxies within a typical dust 
sputtering time tsput ~ 10 yrs, evaluated a few kilopar- 
secs from the galactic center. In addition, since tgput is 
somewhat less than the typical free fall time from sev- 
eral effective radii, if the extended dust were a result of 
merging, ellipticals having far-infrared excesses should 
contain evidence of the merging galaxies within their op- 
tical images, but this is not observed. The mean stel- 
lar ages of elliptical galaxies can be measured from the 
strength of Balmer line absorption features. Annibali et 
al. (2006) find a stellar age of 13.5 ± 3.6 Gyrs for NGC 
4636 but the age for NGC 5044 (14.2 ± 10. Gyrs) is inde- 
terminate because the Balmer absorption is masked by 
the much stronger interstellar Balmer emission in this 
galaxy. The optical colors of both galaxies also seem 
normal, B -V ^ 0.97 - 1.06 (Poulain 1988; Poulain & 
Nieto 1994). Finally, the stellar velocity dispersion and 
rotation are largely normal for both galaxies^ , with 
no kinematical indication of an alien merging stellar sys- 
tem (Caon, Macchetto, & Pastoriza, 2000). We conclude 
there is no optical evidence in either galaxy of a merger 
event during the previous ^ 10* yrs. 

In Temi et al. (2003) we first discovered excess dust 
at 60, 90 and 180 /xm in NGC 4636 from observations 
with the Infrared Space Observatory (ISO) and proposed 
that the short-lived dust resulted from a recent merger. 
However, now that the age of the stars in NGC 4636 
has been found to be very old and the dust is known 
to be extended, our earlier merger explanation is unten- 
able. Instead, as discussed in Temi, Brighenti & Math- 
ews (2007), we now propose that the extended excess 
dust seen in both NGC 4636 and 5044 and other similar 
galaxies has been buoyantly transported from the galac- 
tic core to several kpc following a heating event near the 
central black hole. 

Galactic cores are a natural source of dust. Mathews 
& Brighenti (2003) showed that the - 100 - 300kpc- 
sized nuclear disks that are often observed in the nuclei 
of elliptical galaxies (e.g. Lauer et al. 2005) can result 
from stellar mass loss within the innermost ^^1 kpc at a 
rate Mi « a,Mi « 0.1 Mq yr^^, where Mi ~ 5 x 10^° 
Mq is the stellar mass in the central kpc. 

In the Appendix of Temi et al. (2007) we showed that 
the excess far-infrared emission at 70 and 160 /im cannot 
arise from these small nuclear dust disks/clouds but must 

^ NGC 4636 docs have a few minor pcculicirities. Slit spectra 
of Caon, Macciietto & Pastoriza (2000) indicate that the mean 
stellar radial velocity has a broad central minimum at position 
angle 117°, but not in other directions. Also, for its luminosity 
and central stellar velocity dispersion NCC 4636 has an unusually 
low near-IR stellar surface brightness, causing it to deviate from the 
fundamental plane for galaxy cores (Ravindranath, et al. 2001). 



come from dust distributed 5-10 kpc into the hot gas.^ 
We hypothesized that the spatially extended dust has its 
origin in the nuclear disks/clouds that are disrupted and 
buoyantly transported to large radii in the hot gas as a 
result of intermittent AGN-type heating in the galactic 
core. The mass of excess extended dust, ~ 10^ Mq, is 
consistent with this. 

Consider for example the time tdisk required to form 
an optically obscuring dusty disk of radius r ~ 300 pc. 
To block starlight, the column density in the disk must 
exceed about N ~ 10^^A7'2i cm~^, assuming solar abun- 
dance. So the total mass of an optically dark disk of 
radius r is at least Mdisk ~ [5At/(2 -|- jijlmpNirr'^ « 

3 X 10^Ar2ir3Qo Mq where rsoo is the disk radius in units 
of 300 parsecs and /i = 0.61 is the mean molecular weight. 
[This disk mass is consistent with the upper limit to the 
HI mass in NGC 4636 of 1.1 x 10* M© established by Kr- 
ishna Kumar & Thonnard (1983).] The time for stellar 
mass loss (within the central kpc) to form such an obscur- 
ing disk, tdtsk ~ Mdisk/Mi, is only tdisk ~ 3 x 10'^iV2ir|oo 
yrs, similar to the buoyant rise time in the hot gas. For 
a dust to gas mass ratio 5 ^ 0.01, these disks should 
contain about 3 x 10^ Mq of dust, which is very similar 
to the dust we observe in NGC 5044 and 4636 in excess 
of our steady state model. 

To support this hypothesis further, we show in Figure 

4 Hubble Space Telescope images of the central -^1 — 2 
kpc in NGC 5044 and 4636. Small, chaotically ar- 
ranged dusty fragments are visible against the stellar 
background in both galaxies. The dense gas in these 
fragments is clearly in a highly transient state, dynami- 
cally orbiting in the galactic potential out to a few kpc 
where the freefall time is ~ 10^ years. It is very unlikely 
that these relatively large optically obscuring clouds can 
be interpreted as infalling dust just after being ejected 
from stars. The radius of individual optically dark frag- 
ments in Figure 4 is about 10 parsecs, which translates 
to a dust mass ~ 3 x lO'^ Mq far in excess of that from a 
single star or from a hypothetical agglomeration of sev- 
eral stellar outflows during the short freefall time. We 
propose instead that these dusty regions are fragments 
of a central, more organized disk-like structure that was 
disrupted about lO"" yrs ago by an AGN-related energy 
release. While it may seem unlikely that dense frag- 
ments can result from an AGN energy release and be 
accelerated outwards in this quasi-coherent manner, we 
appeal to studies of Seyfert galaxies where similar large 
and dense clouds are accelerated in just this manner by 
non-thermal, radio-emitting outflows (e.g. Whittle, et 
al. 2005). 

As soon as a dusty disks are destroyed by a significant 
release of AGN energy, another disk begins to form. The 

Cold dust inside dense, optically thick nuclear disks/clouds is 
self-shielded and cannot be heated by the diffuse stellar UV radi- 
ation from hot post-AGB stars. The apparent longevity of these 
optically thick clouds (~ 10^ yrs) suggests that heating by ther- 
mal electrons in the adjacent hot gas is also suppressed. Evidently 
magnetic fields (tangent to the surface of the spinning disks) are 
sufficient to restrict thermal conduction. However, when the dusty 
gas is heated and dispersed by an AGN explosion, some of the dust 
is directly exposed to diffuse galactic UV and hot thermal electrons 
that heat the dust to radiate much more efficiently at infrared wave- 
lengths. But when the same mass of dust is contained in dense, 
self-shielding central disks, its infrared luminosity is unobservably 
low. 
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fragments visible in Figure 4 indicate that some frac- 
tion of the disk received kinetic energy but was not ir- 
reversibly heated by the explosive AGN event. There- 
fore the next generation disks in these galaxies will be 
formed out of these denser remnants of the previous disk 
and new dust that is continuously ejected from red giant 
stars in the galactic core. At present it is impossible to 
know what fraction of disk material is heated to super- 
virial (buoyant) temperatures and what fraction remains 
in the fragments visible in Figure 4. Nor is it clear at 
present what limits the maximum amount of dusty cold 
gas that resides in central disks/clouds. 

The only signs of current nuclear activity in NGC 5044 
and 4636 are small, weak non-thermal GHz nuclear ra- 
dio sources. There is no evidence of disk or Bondi type 
accretion. For example, from Chandra observations of 
NGC 4636 Loewcnstcin et al. (2001) determined that 
the accretion onto the 8 x 10^ Mq central massive black 
hole in this galaxy is less than 3 x 10~^Ledd and less than 
6 X 10""' of the spherical Bondi accretion rate. The gen- 
eral absence of mass accretion in supermassive nucelear 
black holes in elliptical galaxies is consistent with models 
in which most of their mass was acquired during a short- 
lived early evolutionary phase (e.g. Hopkins, Narayan, & 
Hernquist 2006). The AGN- type energy releases in NGC 
5044 and 4636 required in our interpretation must also 
be sufficiently short-lived or radiatively unremarkable. 

4. PAH EMISSION FROM NGC 5044 

Perhaps the most dramatic illustration of extended 
dust in NGC 5044 is the difference image between IRAC 

wavelengths 8 and 4.5 /im shown in Figure 1 and its 
striking similarity to that of the warm gas emission in 
Ha-|-[Nn]. The 8-4.5 /um difference image is constructed 
by subtracting a properly normalized 4.5 /xm image from 
the 8 /Jin image. Emission at 4.5 /xm is due largely to 
stellar photospheres while both photospheric and circum- 
stellar emission contribute at 8 /xm - but the 8 /xm image 
may also contain interstellar emission from hot dust emis- 
sion or the strong 7.7 /um PAH emission feature. When 
the two images are subtracted, both the stellar and cir- 
cumstellar emission are canceled, and only emission from 
the faint interstellar dust emission remains. Our subtrac- 
tion process in NGC 5044 was based on visual examina- 
tion of the surface brightness S in the difference image 
Sg.o— i^S4.5 where the subscripts refer to the IRAC wave- 
lengths and is a variable parameter. As F was varied, 
the stellar surface brightness of the difference image of 
the stellar emission passes through a cancellation when 
F = I]8.o/5^4.5- In this process particular attention was 
paid to the residual stellar brightness in the region along 
the Ha plume about 3-5 kpc south of the galactic center. 
The value of F at which the stellar contributions at 8.0 
and 4.5 microns outside the plume region are equalized 
depends on the galactic radius in NGC 5044. The 8-4.5 
/um image of NGC 5044 shown in Figure 1 corresponds 
to F « 0.52 where the stellar emission near the plume 
canceled and the faint residual interstellar features are re- 
vealed. Our value of F is almost identical to _F = 0.526 
which is the flux ratio F8.0/-F4.5 of a 10 Gyr old solar 
abundance single stellar population estimated by Piovan 
et al. (2003), who include both photospheric and cir- 
cumstellar infrared emission. In our difference image of 
NGC 5044 the interstellar emission exhibits the same ex- 



tended features toward the north and south that are so 
prominent in the Hq!+[NII] isophotes. We estimate that 
the total luminosity of the 8 /um component is about 4 
percent of the total interstellar far-infrared luminosity 
longward of 24 /xm, this is similar to the PAH emission 
fraction observed in Milky Way sources. 

The residual interstellar feature in NGC 5044 could 
be due to the strong PAH emission band at 7.7 /xm or 
to hot dust emission from small, stochastically heated 
dust grains, in some proportion. Dust with temper- 
atures approaching 10^ K could be heated by an in- 
tense beam of optical-UV radiation from the galactic 
center along the direction of the plume. However, to 
our knowledge normal elliptical galaxies with abnormally 
UV-bright nuclei, when viewed along such beams, are 
not commonly observed. In principle, thermal excursions 
of small grains recently heated by individual photons or 
electrons could explain the diffuse 8 /im emission. For 
example, this process is thought to explain the diffuse 
near-infrared dust emission often observed in planetary 
nebulae (e.g. Phillips & Ramos-Larios 2006). During the 
short-lived stages as the dusty hot gas in the southerly 
plume of NGC 5044 cools, the electron density increases 
and stochastic heating of small grains by electron encoun- 
ters may increase. However, in situations where stochas- 
tically heated dust emission dominates, the infrared con- 
tinuum generally increases from 8 toward longer wave- 
lengths since a larger fraction of small grains can be mo- 
mentarily heated to lower temperatures that radiate at 
longer wavelengths (e.g. Draine & Anderson 1985). But 
this increase is not observed in spectra of normal ellipti- 
cal galaxies taken with the Spitzer Infrared Spectrograph 
(IRS) (e.g. Bregman, Temi & Brcgman 2006). Our 24 
fim image of NGC 5044 (not shown in Figure 1) shows 
no evidence of asymmetrically extended emission similar 
to the optical line emission. 

If the extended 8-4.5 feature visible to the south 
in Figure 1 were due to transient thermal excursions of 
small grains, this feature would be even more pronouncieci 
in a 24-4.5 difference image. To explore this possibil- 
ity, we spatially degraded the 4.5 /xm image to the PSF of 
the 24 fim image and varied the flux ratio to cancel the 
stellar and circumstellar emission at both wavelengths. 
We see no non-stellar residual interstellar enhancement 
in the 24-4.5 /xm image of NGC 5044 that is geometrically 
similar to the plume extension visible in optical emission 
and the 8.0-4.5 /i difference image. The absence of ex- 
tended emission in the 24-4.5 /im image of NGC 5044 
is therefore a strong argument for PAH emission at 8 
/xm and for the relative unimportance of thermal excur- 
sions in the mid-infrared dust SED in this galaxy. An 
IRS spectrum of NGC 5044 has been taken but it is not 
currently available in the Spitzer archive, although the 
extended interstellar emission in Figure 1 may be too 
faint to be easily visible as PAH emission features in the 
spectrum. 

In the following discussion we adopt a PAH interpreta- 
tion of the 8 /xm interstellar features in NGC 5044, but 
the precise identification of this dust-related interstellar 
emission is not central to our main arguments. It is un- 
clear if the PAH is located in the warm gas or is in the 
hot interstellar gas where its expected lifetime is short 
but its modes of excitation and emission might be more 
easily excited. 



Extended Dust in Elliptical Galaxies 



5 



Evidence of PAH emission in Spitzer IRS spectra of 
elliptical or other early type galaxies is rare and has 
only been found in a few unusual early-type galaxies that 
have starburst activity in massive cold disks: NGC 4550, 
NGC 4435 (Bressan et al. 2006; both are SBO galaxies) 
and NGC 4697 (Bregman, Bregman & Temi 2006; an E6 
galaxy with a dust lane that may be an SO). 

It is unclear at present if the extended 70 fiia emission 
also shows the same feature to the south of NGC 5044. 
When the 8-4.5 fim feature is convolved with the PRF of 
the 70 /X filter, this southern extension largely disappears. 

The very elongated Hq!+[NII] image in NGC 5044 is 
unusual ~ perhaps the only other elliptical galaxy that is 
known to have a similar feature (unrelated to large disks 
of cold gas) is NGC 5813 (Goudfrooij et al. 1994). In 
normal galaxies without extended dust or warm gas fea- 
tures, all the warm gas is expected to result from stellar 
mass loss and is generally visible only within the cen- 
tral kpc. Gas ejected from evolving stars at a rate a*M* 
(where a, = 4.7 x 10~^° s~^) cannot accumulate in the 
warm gas phase indefinitely since otherwise after ~ 10^ 
yrs the Ha-|-[NII] flux would exceed values commonly ob- 
served in most elliptical galaxies. We conclude therefore 
that the warm line-emitting gas thermally merges with 
the hot (r ~ 10^ K) interstellar gas in this time (Math- 
ews & Brighenti 1999). Therefore the extended warm 
gas observed in NGC 5044 is very transitory, although 
the PAHs and warm gas may be continuously supplied 
for a longer time ^ 10' years as gas cools in the buoyant 
dusty plume. 

5. BUOYANT DUST OUTFLOW IN NGC 5044 

Two questions come to mind about buoyant dust trans- 
port. Can dust grains survive sputtering destruction dur- 
ing their buoyant journey out to ~ 5 kpc, the approxi- 
mate extent of the Ha-I-[NII] and PAH emission features 
in NGC 5044? Can the dust also cool (some of) the buoy- 
ant gas to create the warm gas that emits Ha -I- [Nil]? To 
answer these questions we consider a simple model in 
which gas is assumed to be heated at some small radius 
ro in NGC 5044 and remains in hydrostatic equilibrium 
as it moves buoyantly out at velocity v. For simplic- 
ity we assume v is constant since the buoyant velocity 
depends on the size of the buoyant region and other un- 
certain parameters. After the gas is heated at radius 
ro it has an electron density n that is lower by a factor 
A < 1 than the local ambient density of the hot gas, i.e. 
n(ro) = Ari(ro). In this Section variables with hats refer 
to the (unheated) ambient hot gas in NGC 5044. The 
density and temperature profiles in NGC 5044, found 
from X-ray observations (i.e. Buote et al. 2003), can be 
approximated with the following expressions: 
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and 

G{rkpc) = 1.64(rfepc/30)-°-i°. (4) 

For simplicity we assume that the grains have a single 
initial radius a(ro) and, as they move buoyantly outward, 
are sputtered at a rate given by 
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(Tsai & Mathews 1995) where Up = 0.83ne is the proton 
density. As the buoyant gas flows out in local pressure 

equilibrium, P = P, its temperature follows an adia- 
bat (T oc p('i'~i)/'i' with 7 = 5/3) modified by radiative 
losses, 
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(6) 



The first term in square brackets describes the adiabatic 
cooling of the buoyant gas as it expands, continuously 
adjusting to the decreasing ambient pressure as it moves 
out in the hot gas atmosphere of NGC 5044. The sec- 
ond term is the intrinsic cooling that occurs in the buoy- 
ant gas due to standard radiative losses with coefficient 
^gas{T) (e.g. Sutherland & Dopita 1993) and due to 
the absorption of thermal electrons by inelastic collisions 
with the grains, 
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(7) 

(Mathews & Brighenti 2003). We assume a uniform dust 
mass to gas mass ratio oi 5 = 0.015 in the initial gas, 
similar to that in Milky Way gas. The mass density of 
the grain material Pg = 3.3 g cm^"* is appropriate for 
amorphous silicates. 

The solution to these equations shown in Figure 5 cor- 
responds to the following (somewhat arbitrary) parame- 
ters: ro = 200 pc, A = 0.15, a(ro) = 1/ini. and ?; = 400 
km s^^. The upper panel of Figure 5 shows Chandra 
observations of the hot and cold temperature phases in 
NGC 5044 (Buote et al. 2003). The sohd lines in the up- 
per two panels show the temperature and density profiles 
of the cooler, dense phase T = Tc{r) and the long-dashed 
lines show the buoyant evolution of dust-free gas heated 
by a factor A~-^ at ro- The short-dashed lines in Figure 
5 show the trajectory of the same heated gas but with 
a uniform admixture of dust at a fraction 5 = 0.015 of 
the gas mass. It is seen that the buoyant gas docs in 
fact cool to low temperatures (T -C T) at r « 5 kpc 
similar to the extent of the southern warm gas feature 
observed in NGC 5044. The dust-assisted cooling occurs 
in ^ 10^ yrs even though the buoyant gas was initially 
hotter than the ambient gas T; the bottom panel in Fig- 
ure 5 shows that dust cooling dominates at every radius 
in the dusty buoyant gas. Evidently, the PAH emission 
observed comes from residual dust fragments diminished 
by sputtering. Perhaps the least satisfactory aspect of 
this approximate calculation is the rather high value of 
the buoyant velocity required to reach ~ 5 kpc in ~ 10^ 
yrs, V = 400 km s~-^, which nearly approaches the sound 
speed Cs = {-/kf / ump)^'^ = 480(T/10^ K)V2 km s'^. 
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This suggests that the buoyantly heated gas may have re- 
ceived momentum from a jet emerging from the central 
black hole. 

Considering the many uncertainties involved, it seems 
plausible from Figure 5 that the peculiar asymmetric ex- 
tensions of warm gas in NGC 5044 can be formed by 
dust-assisted cooling of buoyant flow. Buoyant gas that 
was heated less at tq (i.e. A > 0.15) or which con- 
tains a locally higher dust density will cool at smaller 
radius, explaining the full extent of the warm gas feature 
in NGC 5044. The enhanced soft X-ray emission ex- 
tending southward in the general direction of the plume 
(Figure 1) could be due to enhanced emission from cool- 
ing (and therefore denser) hot gas. Indeed, wo note that 
the cooling in Figure 5 occurs on a timescale ~ 10^ yrs 
consistent with other transient events in the hot gas in 
both NGC 5044 and 4636 which we discuss below. 

6. CENTRAL OPTICAL AND X-RAY CAVITY IN NGC 4636 

Optical line emission from warm gas in the central ~ 1 
kpc of NGC 4636 and thermal X-ray emission from the 
hot gas both have a small cavity slightly displaced from 
the galactic center ^ (Figure 6). Demoulin-Ulrich et al. 
(1984), Busonetal. (1994) and Zeilinger et al. (1996), all 
noticed a nearly circular hole of radius rhoie — 1-7" = 120 
parsecs in the Hq:+[NII] image that is clearly visible in 
Figure 6 just NE of the galactic center. At the right 
in Figure 6 is an X-ray image taken from the Chandra 
archives. Superimposed on the X-ray image are approx- 
imate contours showing the radio jets as observed at 1.4 
GHz by Stanger & Warwick (1986) and at 5 GHz by 
Birkinshaw & Davies (1985). The (steep spectrum) dou- 
ble radio jets are about 1 kpc long, oriented at posi- 
tion angle ~ 45°. It is remarkable that the circular hole 
seen by the optical observers is also visible in the Chan- 
dra image. While the X-ray emitting gas (paradoxically) 
appears to have been disturbed very little by the radio 
lobes, the energy source that created the circular cavity 
to the NE has does not contain a significant non-thermal 
radio source of its own. 

Assuming that the circular feature is formed by a point 
explosion, the energy required to form a hole of radius 
Thole in the hot gas is £ « {5P/2p) ■ pVhoie = STrr^^^^P, 
where P = 1.9nkT is the local gas pressure. Using cen- 
tral values of the hot gas density n and temperature T 
from Allen et al. (2006), we estimate that the local gas 
pressure near the circular cavity is P « 6.14 x 10^^" 
dynes cm~^. The energy required to form this hole is 
therefore E f« 200 x 10^^ ergs. This energy far ex- 
ceeds that of a single siipernova, so the cavity was prob- 
ably formed with energy released near the central black 
hole. Since the chaotic warm gas velocities observed by 
Caon, Macchetto & Pastoriza (2000) near the hole are 
V ~ ±50 — 100 km s~^, the age of the hole cannot exceed 
thole ~ '^rhoie/v « 2 X 10'^/?;ioo yrs where vwq is the lo- 
cal warm gas velocity in units of 100 km s~^. Evidently, 
the central AGN-black hole in NGC 4636 has been active 
during the last few million years. 

7. TRANSIENT ACTIVITY IN THE HOT GAS 

The cavity is less visible in the emission line image shown in 
Figure 2 probably because of the gray-scale intensities adopted by 
Caon et aL (2000). 



Unlike most other similar galaxy/groups, XMM X-ray 
spectra of the hot gas in the NGC 5044 group show that 
the hot gas within about 30 kpc from the central galaxy 
does not have a uniform temperature at each radius (Fig- 
ure 5; Buote et al. 2003). In most galaxy groups and 
clusters the hot gas temperature rises from the center 
and peaks at about 0.1 of the virial radius, then decreases 
further out. The emission- weighted (single) tempera- 
ture profile in NGC 5044 has a similar thermal structure. 
However, XMM spectra reveal that NGC 5044 contains 
an additional component of hotter gas in which the tem- 
perature decreases outward from the center (Buote et al. 
2003) as shown by the long-dashed line in the top panel 
of Figure 5. Since this hotter (and therefore less dense) 
gas in NGC 5044 must be buoyant, it is natural to asso- 
ciate it with the dust-transporting regions discussed in 
the previous section. The "multi-phase" character of the 
X-ray spectrum of NGC 5044 is also evident from the 
hot gas density irregularities visible in the Chandra im- 
age (Fig. 1 of Buote et al. 2003) with isophotes shown in 
our Figure 1. It appears that the X-ray feature toward 
the south in NGC 5044 is roughly coincident with the 
Ha-|-[NII] and PAH extensions in that same region. The 
enhanced X-ray emission in this region may result from 
slightly denser gas that is undergoing dust-assisted cool- 
ing in the buoyant plume. Overall, observations indicate 
many dusty plumes in NGC 5044, but only the plume 
toward the south is currently producing warm gas out to 
^ 5 kpc. 

X-ray spectra and images of NGC 4636 also reveal evi- 
dence of recent, possibly somewhat more violent activity 
in this galaxy/group. Chandra observations presented by 
Jones et al. (2002) show a symmetric X-shaped parabolic 
pattern of enhanced emission, presumably caused by 
shock waves, that enclose cavities several kpc to the east 
and west of the galactic center. Jones el al. interpret this 
as the result of two symmetric and nearly simultaneous 
point explosions (~ 10^"^ ergs) offset from the galactic 
center that occurred about 3 x 10^ yrs ago. No radio 
emission has been observed in these cavities which are 
somewhat misaligned from the current radio jets shown 
in Figure 6 (Birkinshaw & Davies 1985; Stanger & War- 
wick 1986). 

The X-ray observations of NGC 4636 and their ex- 
planation do not necessarily require a (buoyant) mass 
outflow that could create the extended dust emission ob- 
served at 70/im. However, strong evidence for outflows in 
NGC 4636 is provided by the XMM spectra and images 
shown by O'SuUivan. Vrtilek and Kempner (2005). Their 
abundance map of NGC 4636 shows a patch of gas with 
high abundance about 25-30 kpc toward the southwest 
of the galactic core. O'Sullivan et al. interpret this as 
evidence of an older (coUimated, plume-like) gas outflow 
that transported gas from the galactic center that was 
enriched by Type la supernovae. Consequently, there is 
now good X-ray evidence that the same kind of transient 
coUimated outflow has occurred in both of these galaxies 
having extended far-infrared emission. 

8. FINAL REMARKS 

We discuss two optically normal elliptical galaxies, 
NGC 5044 and NGC 4636, in which we discovered spa- 
tially extended interstellar dust at 70/xm. In NGC 5044 
we observe asymmetrically extended interstellar S/xm 
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PAH emission associated with warm gas and extending 
about 5 kpc from the center of the galaxy. In view of the 
short sputtering lifetime for extended dust, ~ 10'' yrs, 
and the commonplace optical appearance of both galax- 
ies, we conclude that the extended dust observed cannot 
result from a recent merger with a gas-rich galaxy. The 
stellar morphologies, spectra, colors and old ages of these 
two galaxies appear to be normal with no evidence of an 
additional merger component. Instead, we propose that 
the extended excess dust is created by stellar mass loss 
in the central ~ 1 kpc and stored in small disks or clouds 
in the galactic cores. These disks are intermittently dis- 
rupted and heated by energy released by accretion onto 
the central black holes (AGN feedback). Finally, the 
heated dusty gas is buoyantly transported out to sev- 
eral kpc where it is observed. The physical extent of 
the extended dust is governed by the buoyant velocity 
and either the grain sputtering time or the dust-assisted 
cooling time, whichever is shorter. 

X-ray observations of both NGC 5044 and 4636 show 
unusual, transient features consistent with recent energy 
releases in or near the galactic cores. Both galaxies show 
disorganized, highly fragmented optically visible dust 
clouds in the central 100-200 parsecs that are also highly 
transient. Finally, the timescales for all these processes 
are comparable to the sputtering lifetime of the extended 
dust, about 10^ yrs. 

We are particularly privileged to view NGC 5044 dur- 
ing a brief interval when asymmetric, radially extended 



Allen, S. W., Dunn, R. H., Fabian, A. C, Taylor, G. B. & Reynolds, 

C. S. 2006, MNRAS, 272, 21 
Annibali, F., Bressan, A., Rampazzo, R., Zeilingcr, W. W. & 

Danese, L. 2006, A&A (in press) ( astro-ph/0609175 ) 
Arimoto, N., Matsushita, K., Isliimaru, Y., Ohashi, T. & Renzini, 

A. 1997, ApJ, 477, 128 
Birkinshaw, M. & Davies, R. 1985, ApJ, 291, 32 
Bregman, J. N., Temi, P. & Bregman, J. D. 2006, ApJ, 647, 265 
Bregman, J. D., Bregman, J. N. & Temi, P. 2006, to appear in The 

Spitzer Science Center 2005 Conference: Infrared Diagnostics 

of Galaxy Evolution, ( astro-ph/0604369| 
Bressan, A. et al. 2006, ApJ, 639, L55 

Buote, D. A., Lewis, A. D., Brighenti, F. & Mathews, W. G. 2003, 

ApJ, 594, 741 
Buson, L. M. et al. 1993, A&A, 280, 409 

Caon, N., Macchetto, D. & Pastoriza, M. 2000, ApJS, 127, 39 
Demoulin-Ulrich, M. H., Butcher, H. R. & Boksenberg, A. 1984, 

ApJ, 285, 527 
Draine, B. T. & Anderson, N. 1985, ApJ, 292, 494 
Fazio, G. G. et al. 2004, ApJS, 154, 10 

Goudfrooij, P., Hansen, L, Jorgensen, H. E., & Norgaard-Nielsen, 

H. U. 1994, A&AS, 105, 341 
Hopkins, P. F., Narayan, R., Hernquist, L., 2006, ApJ, 643, 641 
Jones, C., Forman, W., Vikhlinin, A., Markevitch, M., David, L., 

Warmflash, A., Murray, S., & Nulsen, P. E. J. 202, ApJ, 567, 

L115 



8/im PAH and Ha-|-[NII] emission is present in the out- 
flowing gas. We explain this remarkable observation as a 
result of dust-assisted cooling in a buoyant plume of hot 
dusty gas. We show with a simple calculation that dust 
can cool buoyant gas to ^ 10"* K which emits the optical 
emission lines observed. The warm gas phase is main- 
tained in thermal equilibrium near ~ 10^ K by radiative 
losses and UV heating from post-AGB stars. Based on 
the collective Balmer line luminosities of normal ellipti- 
cal galaxies, we conclude that this warm gas can only last 
for ~ 10^ yrs before it thermally melts into the dominant 
hot interstellar phase. Finally, the buoyant outflow de- 
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Fig. 1. — Comparison of optical, infrared and X-ray emission from NGC 5044. The physical scale is 160 parsecs///. The top three panels 
from left to right show Spitzer images at (1) 160 /im, (2) 70 ^tm, and (3) 8 fim. A small source about 20" to the East of the center of NGC 
5044, is revealed in the 8 ^m image. This source is coincident with a point-like source in the B — I image of Goudfrooij et al. (1994) and 
appears as a background galaxy in the HST image. The second row of panels from left to right shows: (1) the difference image 8-4.5 fira, 
(2) isophotes of H«-|-[NII] emission from warm gas in NGC 5044 taken from Goudfrooij et al. (1994) and (3) isophotes from the Chandra 
X-ray images superimposed on an optical image from the Digital Sky Survey. Surface brightness values in the color bars are presented in 
units of ^Jy / arcsec^ . 
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Fig. 2. — Comparison of optical, infrared and X-ray emission from NGC 4636. The physical scale is 73 parsecs///. The top three panels 
from left to right show Spitzer images at (1) 160 fim, (2) 70 fim, and (3) the difference image 8-4.5 ^m. Surface brightness values in the 
color bars are presented in units of fijy /arcsec^ . The second row of panels from left to right shows: (1) isophotes of Ha-|-[NII] emission 
from warm gas in NGC 4636 taken from Caon et al. (2000) and (2) isophotes from the Chandra X-ray images superimposed on an optical 
image from the Digital Sky Survey. 
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Fic;. 3.--- Solid circles show Spitzer MIPS flux observations for NGC 5044 and 4638 at 24, 70 and 160 /im. The lower two thin curves 
show the far-infrared steady state dust SED expected from normal stellar mass loss balanced by sputtering destruction in the hot gas. The 
short-dashed and solid lines show the SED based on initial grain size distributions with maxima at amax = 1 and 0.3 fim respectively. 
These models are clearly insufhcient to explain the observations of either galaxy - an additional source of extended dust is required. Each 
panel shows several models of the dust SED with extra dust, designed to fit the interstellar 70 and 160 /im fluxes (but not 24 /im which 
is circumstellar). As explained in the text, these models are characterized by two parameters {fa, Tex) which have the following values: 
Upper panel: (300,5kpc) (dotted line), (165,10kpc) (long-dashed line), (800,2kpc) [dash-dotted line); Lower panel: (60,4kpc) (dotted line), 
(26,10kpc) (long-dashed line), (133,2kpc) (dash-dotted line). 
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Fig. 4. — Hubble Space Telescope images of the central ~ 1 — 2 kpc in NGC 5044 (left panel) and 4636 (right panel). The optical dust 
maps were generated from archival HST data recorded with the Wide Field and Planetary Camera 2 (WFPC2) using the broad-band filters 
F791W and F547M. We first modeled the stellar contribution, recorded in the broad-band image, with elliptical isophotcs using a two pass 
procedure to properly mask out and remove bad pixels and background/foreground objects. The dust map was then generated from the 
ratio of the continuum image over its purely stellar isophotal model. Dust absorption features axe evident as dark filamentary structures 
in the maps. 
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Fig. 5. — Approximate radial evolution of dust cooling in the constant-velocity buoyant outflow in NGC 5044. Top panel: temperature 
(K) of a dust-free buoyant element (long-dashed line) and dust-filled buoyant element (short-dashed line) compared to that of the ambient 
gas (solid line), Second panel: electron density (cm~^) of the buoyant element (long-dashed line) compared to that of the ambient gas 
(solid line) and that in dusty buoyant gas (short-dashed line), Third panel: time t7 in 10^ yrs for the buoyant element to reach each 
radius (short-dashed line), radius of the grains in microns (short-dashed line). Bottom panel: radiative cooling coefficient Agoa (lower 
short-dashed line) and dust-assisted cooling coefficient A^uat (upper short-dashed line), both in erg cm^ s~^. 
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Fig. 6. — Left Panel: An Ha+[NII] image of NGC 4636 from Zeilinger et al. (1996) also showing the orientation of their slit spectra. 
These authors note that the slits are centered about 1.5" from the galactic center. Right Panel: Chandra X-ray image from the Chandra 
archives with approximate radio contours of the double jet based on Birkinshaw & Davies (1985) and Stanger &; Warwick (1986). The 
contours are symmetric about the green cross that shows the position of the galactic center based on a compromise between images from 
the Digital Sky Survey and the centroid of the HST stellar image. 



